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Abstract. We study single silicon vacancy (SiV) centres in chemical vapour 
deposition (CVD) nanodiamonds on iridium as well as an ensemble of SiV centres 
in a high quality, low stress CVD diamond film by using temperature dependent 
luminescence spectroscopy in the temperature range 5-295 K. We investigate in 
detail the temperature dependent fine structure of the zero-phonon-line (ZPL) 
of the SiV centres. The ZPL transition is affected by inhomogeneous as well 
as temperature dependent homogeneous broadening and blue shifts by about 
20 cm -1 upon cooling from room temperature to 5 K. We employ excitation 
power dependent </' 2 ) measurements to explore the temperature dependent 
internal population dynamics of single SiV centres and infer almost temperature 
independent dynamics. 
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1. Introduction 

In recent years, single colour centres in diamond have attracted wide research interest 
due to exciting potential applications such as solid state based single photon sources 
(for a review see [1]) or solid state spin qubits [2] that can be optically manipulated 
and read out [3j [4J [5]. Silicon vacancy (SiV) colour centres in diamond have been 
identified as a very promising system for the realization of bright, narrow bandwidth, 
room temperature single photon sources: The fluorescence of single SiV centres is 
predominantly concentrated in the purely electronic transition, the zero-phonon-line 
(ZPL), at approximately 740 nm which features a room temperature width of down 
to 0.7 nm [6]. Single SiV centres can be feasibly created by ion implantation [7] or 
by in situ doping during chemical vapour deposition (CVD) [6l[8j|9]. Moreover, the 
emission of single SiV centres is fully linearly polarized at room temperature [5] , which 
is advantageous for applications in quantum cryptography [10] or in the quantum 
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frequency conversion of single photons [IT]. However, the spatial structure of the SiV 
complex as well as its symmetry and the electronic level scheme of the centre are 
still under debate [HJ [T3] . In previous work [T31 [TS] , the SiV centre's ZPL has been 
observed to split into two doublets at low temperature. This unique four line fine 
structure pattern may be regarded as a 'spectral fingerprint' of the SiV centre in high 
quality, low stress diamond and may also be used to unambiguously identify single 
colour centres as SiV centres. The origin of the fine structure is not clear at present; 
in an empirical model it is attributed to a level scheme where excited and ground state 
are split into two sub-states Q3] ■ Depending on the assumed spatial symmetry of the 
SiV complex, the ZPL fine structure has been attributed to a splitting of degenerate 
electronic eigenstates due to a Jahn- Teller effect or to a tunneling between different 
spatial realizations of the centre resulting in the splitting of states [THEI]. Moreover, 
as the spin state of the centre is not clear at present [TB] , one might also expect level 
splittings due to spin-orbit or spin-spin interactions. Due to the unidentified level 
scheme, the use of SiV centres in quantum information applications requiring detailed 
knowledge of the level scheme at low temperature as well as polarization of the optical 
transitions and their linewidth (e.g., [31 [5]) is precluded so far. No detailed studies on 
low temperature or even temperature dependent luminescence of single SiV centres 
are reported in the literature. In this work, we use temperature dependent (5-295 K) 
investigations of single in situ created SiV centres as well as SiV centre ensembles to 
elucidate the ZPL fine structure, line shift and widths. We furthermore investigate 
the temperature dependent population dynamics of single SiV centres. 

2. Experimental 

2.1. Samples 

In the present work, we investigate two types of samples. First, to study the fine 
structure of the ZPL of an SiV ensemble, we employ a single crystal diamond film 
grown using a modified hot filament CVD technique (sample 1). Homoepitaxial growth 
on a (001) lb high pressure high temperature (HPHT) diamond substrate (Sumitomo) 
eliminates stress arising from thermal expansion mismatch of substrate and diamond 
and minimizes the contribution of defects like dislocations. Stress is supposed to be 
the main source of lineshifts and thus inhomogeneous broadening of the SiV centre 
ZPL [17]. We furthermore accomplish high crystalline quality, as witnessed by the 
well defined luminescence properties of the incorporated SiV centres, by applying 
optimized growth conditions involving a low methane fraction (0.26% CH4 in H2) 
and slow growth. Using the same parameters for the growth of an additional sample 
on a non-diamond substrate, we estimate a growth rate of approximately 10 nm/h 
equivalent to a total film thickness of 80-100 nm. SiV centres are created in situ due 
to residual silicon contamination of the CVD reactor. We point out that the HPHT 
substrate does not contain SiV centres, thus the observed SiV centre ensemble fully 
resides in the high quality CVD diamond film. 

Second, we employ spatially isolated nanodiamonds (NDs) grown using microwave 
plasma CVD on an iridium (Ir) surface to investigate single SiV centres. The 150 
nm thick Ir layer is deposited onto a thin buffer layer of yttria-stabilized zirconia 
on a silicon wafer [18] . Randomly oriented NDs are grown on seed diamonds 
(Microdiamant, MSY GAF 0-0.03 /im) which are distributed on the substrate by 
spin coating from an aqueous solution (for details of the seeding process see [6]). We 
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use two ND samples: Sample 2a contains NDs with an average size of 130 nm (growth 
conditions: 25 min, 1% CH 4 in H 2 , 2000 W microwave power, 30 mbar); Sample 2b 
contains NDs with an average size of 220 nm (growth conditions: 55 min, 0.4% CH4 in 
H2, 2000 W microwave power, 30 mbar). Scanning electron microscope (SEM) images 
of sample 2a can be found in |6]; figure Q] displays SEM images taken from sample 2b. 
By adjusting the seed density, we obtain approximately 2 NDs per /1m 2 on sample 2a 
and approximately 0.3 NDs per /j,m 2 on sample 2b. In both samples, single SiV centres 
are created in situ due to plasma etching of the underlying silicon substrate at its 
uncovered edges and subsequent silicon incorporation. As discernible from figure[T]and 
reference [6], both samples display highly facetted NDs featuring shapes close to the 
cubo-octahedral shapes of single crystal diamonds. The growth conditions employed 
for sample 2b seem to be superior compared to the conditions used for sample 2a as 
in the SEM images the NDs are found to exhibit even cleaner surfaces without cracks 
and distortions. In contrast to the homoepitaxial film, individual NDs, in general, 
exhibit a variable stress state that develops during the growth (for details see, e.g., 
[6l [9]). The NDs on Ir have been grown using methane with a purity of 99.9995% 
and hydrogen of initially 99.999% which was further purified by a silver palladium 
membrane in order to reduce the nitrogen content by several orders of magnitude. 
We also note that no nitrogen vacancy centre luminescence could be detected. The 
nitrogen electron-paramagnetic-resonance signal in bulk crystals formerly grown under 
these conditions was below the detection threshold («40 ppb). 



(a) (b) 




100 nm 



Figure 1. SEM images of CVD grown NDs on sample 2b: (a) overview image 
displaying spatially isolated NDs (b) detailed images of well faceted NDs. 

2.2. Experimental Setup 

To analyze the colour centre fluorescence, we employ a homebuilt confocal laser 
microscope. For the temperature dependent single emitter investigations, a flow 
cryostat (Janis Research, ST-500LN) operated with liquid helium is integrated into the 
optical setup. Temperatures between approximately 5 K and room temperature (295 
K) can be selected via resistive electrical heating of the cryostat sample mount. The 
temperatures given throughout this paper have been measured at the sample mount of 
the cryostat. The excitation laser, either a frequency doubled diode pumped solid state 
laser (at 671 nm or 532 nm) or a tunable Tksapphire laser (here 690 nm), is focused 
through a microscope objective with NA=0.8 (Olympus, magnification lOOx). Using 
a dichroic beam splitter and dielectric filters, we separate the SiV fluorescence in the 
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spectral range 730-750 nm from residual laser light reflected from the sample surface. 
The fluorescence is coupled into a multimode fibre. Subsequently, the fluorescence 
is directed to a Hanbury Brown Twiss (HBT) setup (avalanche photodiodes, APDs: 
Perkin Elmer SPCM-AQRH-14, correlation electronics: Pico Quant, Pico Harp 300) to 
measure the intensity autocorrelation (g( 2 ) ) function of the fluorescence. Alternatively, 
the fluorescence can be analyzed spectrally. We first use a grating spectrometer 
(Horiba Jobin Yvon, iHR 550). For cryogenic temperature spectroscopy, we employ 
a grating with 1800 grooves/mm enabling an approximate resolution of 0.05 nm (30 
GHz) at 738 nm. 

To perform emission spectroscopy with higher spectral resolution the setup is 
also equipped with a scanning Fabry-Perot interferometer (FPI). The FPI exhibits a 
free spectral range (FSR) of 34 GHz in a planar mirror configuration and a hnesse 
F of 54 (yielding a resolution of 0.6 GHz). The comparably low finesse was chosen 
intentionally in favour of a high signal transmission of T = 86%. We preselect a single 
line out of the four line fine structure using a silica etalon with dielectric coating 
(reflectivity R = 94%) which acts as a 44 GHz bandpass filter. The preselection is 
necessary as the fine structure of the ZPL spans a spectral range much larger than 
the FSR of the FPI. The transmitted fluorescence is detected using the APDs of the 
HBT setup. 

3. Results 

3.1. Photoluminescence spectra 

First, we investigate the low temperature photoluminescence spectra of ensembles 
of SiV centres. As a preliminary remark, we note that in polycrystalline diamond 
films on silicon (sample description see, e.g., [8]), we recently found a temperature 
independent, inhomogeneous ZPL width of approximately 2 nm that easily dominates 
a ZPL fine structure. This observation is in accordance with the literature [19j [14] 
and illustrates the sensitivity of the SiV ZPL wavelength to the local environment. 
In contrast, homoepitaxial CVD diamond films as introduced in section |2~T1 (sample 
1) are a host material with high crystal quality and low stress level enabling the 
observation of the fine structure of the ZPL of an SiV ensemble as displayed in 
figure Ela): At temperatures below 70 K, the four line fine structure develops, 
accompanied by additional peaks stemming from less abundant silicon isotopes (Si 29 , 
Si 30 ) in accordance with the literature [141 ITS] . This fine structure has been explained 
using an empirical level scheme with a split ground and excited state [14] . We 
determine the (larger) excited state splitting to be 260 GHz, whereas the ground 
state splitting amounts to 50 GHz, thus coinciding with observations in the literature 
within 10%. Upon cooling, the ZPL blue shifts by approximately 1 nm (18.4 cm -1 , 
2.3 meV) reproducing the observation in [TS] within 20%. The observed linewidth 
of the luminescence spectrum at low temperature is limited by the resolution of 
the employed grating spectrometer. To further investigate the linewidth, we thus 
employ high resolution emission spectroscopy using a scanning FPI. The obtained 
data are shown in figure [2jb) for each of the four fine structure lines. The measured 
linewidths range between 9 GHz and 15 GHz. Similar values have been reported for 
the inhomogeneously broadened ZPL of an ensemble of NV centres in high-quality 
CVD-grown diamond [20]. As the inhomogeneous distribution of ZPL wavelengths is, 
in general, larger for the SiV centre than for the NV centre (see e.g. [9]) the linewidth 
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Figure 2. (a) Temperature dependence of the ZPL fluorescence in the 
homoepitaxial CVD diamond film (sample 1) excited with 532 nm laser light. 
Spectra up to 150 K allow for the observation of the temperature dependent 
variations of the fine structure. The inset shows the ZPL at room temperature for 
comparison, note the different wavelength scales, (b) High resolution spectrum 
of the fine structure lines at 4.8 K recorded using a scanning FPI (black dots, 
excitation 690 nm). Note that the spectrum shown here is composed of the 
measurements of the individual fine structure lines preselected as discussed 
in section 12.21 The curves have been fitted using Gaussian peak functions 
(red solid lines) to yield the optimal correspondence with the measured data. 
The underlying grey spectrum, acquired with the grating spectrometer, serves 
as a guide to the eye to illustrate the correspondence between the FPI and 
spectrometer measurements. 



measurements obtained here indicate the high structural quality of the homoepitaxial 
diamond layer. 

In spatially isolated CVD nanodiamonds on Ir, investigation of single, bright 
SiV centres is feasible [5JIH]- We here discuss four emitters in detail, labeled C1-C4. 
Emitter CI is hosted by a ND on sample 2a, the other emitters are found in sample 2b. 
Figure [3] summarizes temperature dependent luminescence spectra of the investigated 
emitters. Emitter C1-C3 show a four line fine structure, partially with badly resolved 
components, while emitter C4 displays a spectrum dominated by a single line. Prior 
to discussing the spectra in detail, we address the question wether these multi-line 
spectra, considering the significantly varying brightness of single SiV centres [22 , can 
be clearly identified as single emitter spectra. Section [3T4l discusses g( 2 > measurements 
of the SiV fluorescence under continuous laser excitation. We here anticipate some 
of the results: the g^ measurements indicate the presence of single emitters; the 
measured deviations from g( 2 '(0) = are almost fully explained by the timing jitter 
of the measurement setup and in some cases luminescence background. However, we 
find residual deviations < 0.15. These further deviations from the ideal case possibly 
originate from the presence of a second weak emitter. Assuming that a residual of 
</ 2 )(0) of 0.15 is unambiguously detectable, we calculate that it is feasible to identify 
emitters with an order of magnitude difference in their brightness. Examining the 
relative line intensities, g^ measurements for emitter CI and C3 cannot exclude that 
the weakest line component originates from a second emitter. For emitter C2, the 
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Figure 3. Temperature dependent photoluminescence spectra of single SiV 
centres under 671 nm laser excitation: (a) emitter CI (0.32 mW excitation power), 
(b) C2 (0.16 mW), (c) C3 (0.06 mW below 30 K, 0.6 mW above 30 K) and (d) 
C4 (0.1 mW below 10 K, 0.07 mW above 10 K). 



same holds for the shorter wavelength doublet. However, the following analysis of the 
fine structure pattern also aids in assigning the observed lines to single SiV centres 
and their level schemes. 

The observed single emitter spectra display three major differences compared to 
the spectrum of the SiV ensemble in the homoepitaxial CVD diamond film: First, 
the absolute positions as well as the spacings and the relative intensity of the fine 
structure components are altered and vary between different emitters. Second, the 
room temperature linewidths are significantly smaller for the single emitters. Third, 
the lines at low temperature display additional broadening. The change of the room 
temperature linewidth will be addressed in section [3.2.1 1 the additional line broadening 
will be discussed in section 15.2.21 The varying spectra arise from changes in the local 
environment of the individual colour centres: Stress in the diamond lattice generally 
changes the energy levels of colour centres. Symmetry conserving stress changes 
the absolute energetic position of the levels without affecting the level splittings, 
whereas symmetry breaking stress changes the level splitting [23]. For SiV centres, the 
splittings AE e , AE g of the excited and ground state, respectively, have been shown 
to vary under externally applied stress in (f 00) direction: for a stress of 0.4 GPa, AE e 
amounts to 4.5 times AE e without stress, whereas AE g even increases to a value of 11 
times AE g without stress [T3]. These changes lead to altered positions and spacings 
of fine structure components and are accompanied by relative intensity changes of the 
fine structure components [15]. Particularly, already under a moderate stress of 0.33 
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GPa in (100) direction, a single line dominates the spectrum while other fine structure 
components mostly vanish [TQ. Thus, even the strongly modified spectrum of emitter 
C4 can be attributed to an SiV centre with a ZPL fine structure altered by local stress 
fields. We note that we also observe emitters where only a doublet line is present in 
the spectrum. For emitter C2 (CI), the fine structure pattern is comparable to the 
pattern of the ensemble of SiV centres in homoepitaxial diamond: Despite the fact 
that AE g is twice (three times) the spacing for the ensemble, the ratio AE e /AE g 
matches the ensemble case within 6% (22%). The results in [15] document a change of 
AE e / AE g for (100) stress. However, no information on the response of SiV centres to 
stress in other directions is found in the literature. Furthermore, the stress distribution 
inside the investigated NDs is unknown. We thus attribute the variety of observed 
line patterns, for which the discussed emitters are just examples, to the local stress 
states experienced by the individual colour centres. Moreover, static electric fields due 
to charges in the crystal might influence the spectrum. Data on the Stark shift of the 
SiV centre is not available in the literature. 




735.5 736.0 736.5 737.0 737.5 730 731 732 733 734 735 736 
Wavelength (nm) Wavelength (nm) 



Figure 4. Low temperature spectra for different emitters demonstrating 
instabilities [(a) 8 K, (b) 30 K]. All spectra are normalized to the peak value 
and displaced for clarity, (a) spectrum of an emitter with resolution limited lines, 
but unstable relative intensity (integration time 60 s blue and red curves, 180 
s black curve), (b) emitter with spectral jumps, but similar spectral pattern in 
both states (integration time 5 s). 

Whereas many of the SiV centres show stable emission, some emitters with 
unstable emission spectrum were observed. Figure HJa) displays an emitter showing 
narrow, spectrometer resolution limited lines. Remarkably, however, the relative 
intensity of the fine structure lines varies over time. Figure [4jb) , in contrast, displays 
spectra taken from an emitter undergoing spectral jumps. Each spectrum has been 
integrated for 5 s and witnesses spectral jumps on a timescale of seconds: Spectra 
revealing one (both) of the ZPLs arise if the emitter resides in one (both) of the 
states during the spectrum integration time. Remarkably, the spectra in both states 
resemble a very similar ZPL fine structure despite the significant wavelength difference. 
The jumps might occur due to the transfer of an electron between SiV centres under 
(significantly) different stress conditions if luminescence is only observed from the 
momentarily charged centre. Spectral jumps have also been reported in 24 and 
[25] for colour centres emitting around 770 nm and chromium-related centres. In 
accordance with [2TJ [26l [24] , we furthermore observe blinking or bleaching of single 
centres. The emission instabilities as well as the additional broadening, commonly 
referred to as spectral diffusion, might be connected to charged defects or fluctuating 
charges inside or at the surfaces of the NDs: 
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Ionization of other impurities inside the NDs A prevalent impurity in diamond is 
substitutional nitrogen [N s ] which forms a deep electron donor with an ionization 
energy of approximately 1.7 eV (729.3 nm, e.g., [27]). Thus, nitrogen donors in the 
focal volume of our 671 nm or 532 nm excitation laser might be photoionized, forming 
[N s ] + and electrons in the conduction band. Additionally, very recent experiments 
show that [N s ] can also act as an acceptor [28]. Due to these processes, fluctuating 
charges in the conduction band of diamond are created and shift colour centre 
transitions via the Stark effect [27l [29] . For single NV centres, linewidths of up to 90 
GHz in nitrogen rich type lb diamond have been observed [37] • As the [N s ] content 
of the NDs has been estimated to be below ps 40 ppb (see section [2~T]) . wc assume 
that only a few nitrogen atoms - if any - are incorporated in NDs with a size of 130 
nm. For the homoepitaxial film, we cannot estimate the nitrogen content. However, 
it is supposed to be significantly higher than for the CVD NDs. It should be noted 
that the lb HPHT diamond substrate used for the growth of the homoepitaxial film 
as well as the HPHT seed diamonds used for the ND growth contain a significantly 
higher concentration of [N s ] (in the order of 100 ppm, [17]). As a further effect, 
other impurities in diamond exist in multiple charge states and switch their charge 
state continuously under laser excitation (photochromism) . Charge state conversion 
has been observed for NV centres [30^ as well as for SiV centres [16]. Thus, also 
the charging or discharging of other colour centres may contribute to a fluctuating 
charge background and might lead to line broadening or unstable emission. It remains 
unclear, why the ensemble in the presumably less pure homoepitaxial diamond is less 
prone to line broadening and emission instabilities. 

Trapping and release of charges on the surfaces of NDs Recent work on 5 nm 
detonation NDs identifies surface charge traps due to impurities, dislocations and 
dangling bonds as the major cause of fluorescence blinking of NV centres |31| . 
Similarly, NV centres implanted close to the surface of CVD diamond show strong 
spectral diffusion compared to in situ formed centres deeper in the same sample |32) . 
These findings suggest that also light induced fluctuating charges [33] at surfaces 
induce spectral diffusion and non-stable emission. However, we point out that the SiV 
centres in the 80-100 nm thick homoepitaxial film also experience a close proximity 
to the surface without suffering from spectral diffusion. 

Influence of the Ir substrate In addition to free charges created in the diamond, 
one might think of charge transfer between the diamond and the metal substrate. 
Ir(001) has a workfunction of 5.67 eV [34]. Assuming a hydrogen terminated diamond 
surface (band diagram see [35]), electrons from the valence band in diamond might be 
transferred to unoccupied states of the Ir metal by visible light (density of states for 
Ir see [36]) and thermalize back to the valence band of diamond. Thus, also the light 
induced charge transfer between substrate and diamond might add to a fluctuating 
charge background und thus to spectral diffusion of SiV centres. 

3.2. Temperature dependent emission linewidth 

The emission linewidth of single colour centres is a crucial parameter considering 
their application in quantum information processing, in particular considering the 
generation of indistinguishable photons. The spectra of figure [3] reveal a temperature 
dependent broadening of the ZPL. In general, both inhomogeneous and homogeneous 
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line broadening occur, where the latter relates to dephasing processes due to phonons 
[23] . The temperature dependent homogeneous broadening is a results of quadratic 
electron phonon coupling and leads to a Lorentzian lineshape Jl(A) [37] [35] • Especially 
for single emitters in NDs, spectral diffusion can induce additional broadening. This 
mechanism constitutes an inhomogeneous broadening and results in a Gaussian 
emission line shape /g(A) [21] • If the charge fluctuations discussed above are 
temperature independent, one would assume spectral diffusion to be temperature 
independent, as observed for NV centres in [33] in the temperature range 5-20 K 
and in |39j for chromium related centres. 

If both line broadening types occur the emission line shape is given by the 
convolution of a Gaussian and a Lorentzian line. This line shape is termed a Voigt- 



icg and u>l are the Gaussian as well as Lorentzian full widths at half maximum 
(FWHM) of the emission line, A c is the peak wavelength and A is proportional 
to the area of the line. By fitting a Voigt-profile to the measured line shapes, we 
separate homogeneous and inhomogeneous line broadening contributions for the ZPL. 
The resulting curves are displayed in figure [3j and match the measured data very 
well. Note that for the following discussion we have corrected the obtained linewidths 
for the spectrometer response function (given also by a Voigt profile). A general 
trend considering the line broadening is observed: For the SiV centre ensemble in 
homoepitaxial CVD diamond (sample 1), at 5 K linewidths of below 20 GHz have been 
found (see section I3~T|) . For emitters C1-C4 in contrast, figure [3] reveals significantly 
broadened lines (20-150 GHz) at around 5 K with lineforms close to a Gaussian shape. 
For elevated temperatures, the observed lineform changes and is closer to a Lorentzian 
line. 

Considering the evaluation of the data, the complex spectral pattern together 
with the large number of free parameters for each Voigt-profile renders least-chi 2 
fitting of the emission spectra extremely challenging. Therefore, we have adapted the 
following strategy to fit the line patterns: At low temperature, we assume wl to be 
far below the spectrometer resolution; we thus fix wl to the value obtained from the 
spectrometer lineform. The fits result in a constant value for wg up to about 50 
K. After correction for the Gaussian contribution of the spectrometer response, wq 
gives the spectral diffusion linewidth for the respective line. We assume wq to be 
temperature independent as shown, e.g., in [39 for chromium related centres. We 
thus fix wq to the (mean) value obtained for each line at low temperature and fit wl 
representing the increasing homogeneous broadening for higher temperature. 

3.2.1. Homogeneous linewidth The temperature dependent homogeneous linewidth 
is determined by the phonon spectrum of the perturbed host crystal, i.e., the crystal 
including the colour centre, as well as the population of the phonon modes (for details 
see [37]). The phonon spectrum of the perturbed crystal cannot be directly measured, 
however, the line broadening can be approximately related to the one-phonon assisted 
absorption spectrum of the colour centre [37]. Measuring the latter is not feasible in 
our experiment. Furthermore, deriving the absorption spectrum from the emission 
spectrum is precluded due to a lack of inversion symmetry of the SiV spectrum about 
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the ZPL [IT]. We thus cannot use this approach for a qualitative analysis of the line 
broadening. 

Treating the crystal in Debye approximation, Maradudin derives a T dependence 
of the linewidth for low temperatures if only non-degenerate electronic states take 
part in the transition [37]. For degenerate electronic states, phonon interactions can 
induce transitions between degenerate sublevels of electronic states (dynamic Jahn- 
Teller effect). Thus, in the case of degenerate electronic states, phonon interactions 
can result in additional dephasing mechanisms. Hizhnyakov et al. extend the 
theory to degenerate electronic states and include the effect of a softening of bonds 
(rearrangement of electrons) and thus an enlarged low energy local phonon density 
in the excited state of the colour centre [42]. They find a aT 3 + cT 7 (a,c = const.) 
dependence and verify for an NV centre ensemble a dominant influence of the T 3 
dependence. On the other hand, Hizhnyakov and Reineker [33] attribute a T 3 
dependence to fluctuating fields which are created as a phonon modulates the distance 
between the colour centre and other defects in the crystal. The authors emphasize 
that the T 3 dependence has been observed in systems exhibiting a significant 
inhomogeneous broadening of ensembles, thus indicating an influence of the purity 
of the host material on the homogeneous broadening. Using perturbation theory, a T 5 
dependence of the linewidth at low temperature in the presence of a dynamic Jahn- 
Teller effect has been derived [32] • In principle, taking into account the discussion 
above, it should be possible to identify the main homogeneous line broadening 
mechanism from the dominant term in the polynomial aT 3 +bT 5 +cT 7 (a, b,c — const.) 
dependence of the linewidth on temperature. Nevertheless, observations for NV 
centres in the literature illustrate that the results are not at all unambiguous: Either 
dominant T 3 , T 5 or T 7 broadening has been reported [3H H3J H3 EH- Figure [Sf a) 
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Figure 5. Homogeneous linewidth of (a) the SiV ensemble in a homoepitaxial 
CVD film (sample 1) as well as (b) selected single emitter (C3). The main 
peaks are used for the analysis as they can be observed over a wide range of 
temperature. The peaks are consecutively numbered with increasing wavelength. 
For the ensemble (single emitter) the fits take into the account the data from peak 
2 (peak 3). 

displays the temperature dependent homogeneous linewidth for the ensemble of SiV 
centres in the homoepitaxial CVD diamond film; Figure [5Jb) shows the same for 
emitter C3 as an example for a single emitter. A difference in room temperature 
linewidth is clearly visible. Fitting the temperature dependent linewidth to either a 
T 3 , T 5 or T 7 dependence, it is clear from figure[5]that the latter does not at all describe 
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the experimental data. For the ensemble in the single crystal as well as emitter CI, C3 
and C4, the T 3 dependence gives the lowest % 2 , while only for emitter C2 T 5 leads to 
the best fit. Taking into account the discussion above as well as the observation of a 
significant inhomogeneous line broadening (spectral diffusion) we tentatively attribute 
the main line broadening mechanism as related to defect induced broadening. We note 
that this observation is similar to very recent results obtained using single chromium 
related colour centres [39] . 



3.2.2. Spectral diffusion For the majority of emitters in NDs, we observe additional 
inhomogeneous broadening attributed to spectral diffusion at low temperatures. 
Figure [B] summarizes the Gaussian spectral diffusion linewidths of emitters C1-C4. 
For emitter C3, only the short wavelength peak and the main peak of the red shifted 
doublet are taken into account. Note that badly resolved peaks and an unclear peak 
structure lead to significant errors for this emitter. For emitter C4, we consider only 
the main peak. Emitter C4 noticeably shows the lowest spectral diffusion linewidth, 
despite the strongly modified spectral pattern. For all emitters, the measured 
linewidths vary between 25 GHz and 160 GHz. Significant differences between 
individual lines of the same emitter may occur with linewidth differences up to nearly 
50%, indicating that different electronic transitions can involve different susceptibility 
to spectral diffusion. The difference between individual emitters indicates a strong 
dependence on the local environment, i.e., proximity to possibly ionized impurities 
or surfaces as also reported for NV centres in NDs [26]. From the measurements 
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Figure 6. Spectral diffusion linewidths of single SiV centres in NDs. The 
errors reflect the standard deviation of the width obtained at low temperature. 
As discussed in section I3.2I we use low temperature values to estimate the 
spectral diffusion linewidth up to the temperature where a systematic increase 
in homogeneous linewidth is found (CI: f 10 K, C2: 45 K, C3: 80 K, C4: 10 K). 

discussed in section I3~4l we estimate the excited state lifetime of emitters CI, C3, 
C4 to range from 0.72-1.28 ns yielding a lifetime limited linewidth Ai/ft of 0.12- 
0.22 GHz. Thus, for emitter CI and C3, the smallest spectral diffusion linewidth 
corresponds to 450 • Auft and 170 • Avft for emitter C4. We note that we also 
observe spectrometer resolution limited lines of single SiV centres in CVD NDs [see 
figure Ufa)]. Furthermore, we find a linewidth of approximately 10 GHz for an SiV 
ensemble in a high quality, homoepitaxial CVD film (sample 1, see section l3"TTj) . This 
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indicates that narrow lines can be obtained for SiV centres in low-defect environments. 
We suggest that the main cause of spectral diffusion is the non-resonant excitation 
laser as indicated by the following observation and also reported for NV centres in the 
literature, e.g., in [331 113 US] . For emitter C4, the inhomogeneous linewidth changes 
from 20 GHz at 0.3 Psat to 55 GHz at 14 Ps a t- We thus suggest that by using resonant 
excitation, reducing the laser photon energy as well as the necessary excitation power, 
lower linewidths are feasible for single SiV centres in NDs. 

3.3. Temperature dependent ZPL shifts 




Temperature (K) Temperature (K) 



Figure 7. ZPL line shift (a) for the ensemble of SiV centres (b) for single 
SiV centre (emitter C2). The peaks are consecutively numbered with increasing 
wavelength. Note that the fit takes into account the data of peak 3. 

For the application of SiV centres as sources of indistinguishable photons as well 
as for the coupling of SiV centres to resonant photonic structures, knowledge of the 
precise temperature dependent wavelength of the ZPL is crucial. All observed SiV ZPL 
fine structure components show a similar blue shift upon cooling. Figure [7] displays 
the shift for one peak of emitter C2 as well as for the ensemble in homoepitaxial CVD 
diamond. The overall shift ranges from -19 cm -1 to -31 cm -1 . Thus, the observed 
blue shift is in accordance with reports in the literature (-26 cm -1 , [49]; -22 cm -1 , 
[T9]). Note that we choose intense peaks for the analysis as observation of the weak 
fine structure components is precluded at elevated temperatures; furthermore, the 
coalescence of peaks leads to uncertainties at elevated temperatures. 

Blue shifting of the ZPL upon cooling has also been observed for other colour 
centres, e.g., NV centres [42j [45] . The cause of the lineshift as well as the exact 
temperature dependence are discussed controversially in the literature: During cooling 
the diamond lattice contracts (for recent measurements of the temperature dependent 
lattice constant see [SD]). Lattice contraction as well as temperature dependent 
electron lattice interaction lead to a temperature dependent bandgap in diamond [51] . 
Feng et al. [19] relate the ZPL shift to a change in the bandgap Eg of diamond (for 
details see [E]). However, the change of Eq measured for diamond [51] is not directly 
proportional to the ZPL shift observed here. Thus, we conclude that the change of 
Eg cannot be straightforwardly used to describe the temperature dependence of the 
ZPL shift without making further assumptions. 

Lattice contraction furthermore introduces lineshifts similar to shifts due to 
mechanical stress. The ZPL shift AEi induced by lattice contraction can be calculated 
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where c\\ and C12 are the elastic moduli of diamond, a v is the volume expansion 
coefficient of diamond and A is the shift rate of the ZPL under hydrostatic compression. 
As the parameter A is not reported in the literature for the SiV centre, we cannot 
quantify AEi. Reference [50] finds empirically that below approximately 100 K the 
temperature dependence of a v matches a T 3 dependence, leading to a T 4 dependence 
of AEi. Figure [7] shows the T 4 dependence fitted to the experimental data. It 
reasonably well describes the temperature dependent shift for the single emitter but 
shows deviations at temperatures above 200 K whereas the shift for the ensemble is 
only moderately described. Note the larger error of the shift at elevated temperature 
due to the coalescence of line components. 

As discussed in section 13 .2.11 the investigated SiV centres exhibit a temperature 
dependent homogenous linewidth. Thus, a quadratic electron-phonon coupling is 
present [38] that also contributes to the lineshift. The overall lineshift is given as 
the sum of the phonon dependent parts and the lattice contribution AEi. Without 
knowledge of the parameter A, we cannot separate the lineshift due to phonon coupling 
and lattice contraction. The ZPL shift due to quadratic electron-phonon coupling has 
been derived to be proportional also to T 4 for temperatures well below the Debye 
temperature T D of the solid [37] (T D w 1900 K for diamond 42 ). Additionally, in 
analogy to the temperature dependence of the linewidth, one can use the one-phonon 
absorption spectrum together with the phonon population to calculate the line shift. 
Due to the same line of argument as for the line broadening, using this approach is not 
feasible to evaluate our data. Furthermore, Hizhnyakov and coworkers l 42| predict a 
aT 2 + bT 4 (a, b — const.) dependence when taking into account the softening of bonds 
in the excited state (see also discussion in section I3.2.1[) . We also fit a aT 2 + bT 4 
dependence to the measured data which more closely fits the results as compared 
to the T 4 model (see figure [Jj. Nevertheless, as Hizhnyakov and coworkers do not 
comment on the influence of the lattice expansion in their work, this model does not 
allow for the separation of lattice and phonon contribution either. Summarizing, the 
measured ZPL shift is consistent for different emitters and fine structure components 
and the overall shift agrees with previous observations on SiV centres. The measured 
data is consistent with a aT 2 + bT 4 temperature dependence. However, we cannot 
identify the individual contribution of lattice contraction and phonon coupling to the 



3.4- Temperature dependent population dynamics 

We here deduce the full population dynamics of single SiV centres at different 
temperatures from excitation power dependent measurements of the g^ function 
as discussed below. To estimate the temperature dependence of the dynamics, we 
perform measurements at room temperature and cryogenic temperature (20-30 K). 
The g( 2 ^ measurements prove single emitter nature for all observed SiV centres. The 
brightness of single emitters is deduced from saturation measurements. The observed 
maximum photon rates 1^ and saturation powers P sa t are summarized in table [T] (for 
a general discussion of saturation curves see, e.g., [21]). Note that the values at low 
temperature have been corrected for a lower excitation and collection efficiency of the 
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confocal setup due to imaging through the cryostat window. 1^ for emitters CI and 
C3 agrees within less than 10% at room and cryogenic temperature. I x for emitter C4 
changes considerably, indicating a higher count rate at low temperature. However, we 
point out that for emitter C4, background fluorescence had to be taken into account 
and was especially pronounced at room temperature. Thus, over-accounting for the 
background might explain an underestimated value of /oo at room temperature. In 
the literature, temperature dependent non-radiative quenching of SiV centres has been 
reported: Feng et al. [19] report a decrease of SiV luminescence by a factor of 8 between 
50 K and room temperature. Sternschulte et al. [TS] find that the excited state lifetime 
of SiV centres in bulk diamond decreases from 4 ns to 2.7 ns upon heating from 5 
K to 300 K, indicating possible lifetime shortening due to non-radiative quenching. 
For single emitters, an enhanced quantum efficiency due to the suppression of non- 
radiative decay channels might influence the maximum obtainable photon rates at 
low temperature. From our observations, however, we infer that no significant de- 
activation of non-radiative channels occurs. This conclusion is also supported by the 
temperature dependent g( 2 > measurements discussed below. The observed saturation 
powers Psat are of the same order of magnitude for room and cryogenic temperature. 
Thus, we exclude any significant influence of the temperature on the off-resonant 
excitation process. 



Emitter 


Psat RT 


loo RT 




Psat 


loo 


APsat 




CI 


73 (iW 


2.39 Mcps 


25 K 


97 /J.W 


2.17 Mcps 


-33% 


9% 


C3 


65 /xW 


0.775 Mcps 


20 K 


57 fiW 


0.713 Mcps 


13% 


7% 


C4 


32 


0.593 Mcps 


30 K 


18 fiW 


1.08 Mcps 


45% 


-82% 



Table 1. Saturation power P sa t and maximum photon count rate Too for emitters 
CI, C3, C4 at room temperature and at cryogenic temperature. The temperature 
for the cryogenic measurements is given in column Ti ow . AP sa t and Aioo give 
the deviation of the low temperature value from the room temperature value in %, 
calculated from AToo = [loo(RT) — I<x>{Ti ow )]/ Ieo(RT) and analogous for AP SB t. 

Figure [5] exemplarily displays the excitation power dependent </ 2 ' functions for 
emitter CI at room temperature and 30 K. The measured g^ data are fitted with 
the <?( 2 ) function of a three level system 

3 (2) (r) = 1 - (l + a)e-l T l/ Tl +ae-l T l/ T2 (3) 

which we furthermore modify to take into account the instrument response and 
possible background fluorescence (for details see The employed three level 

system, described in detail in [SJ [3T], takes into account two levels for the ZPL 
transition (1 and 2) and a shelving state (3). The fits agree well with the measured 
data for all delay times r and identify all SiV centres as single emitters. The deviation 
of <?( 2 )(0) from is induced by the instrument response of the measurement setup as 
confirmed by the accordance of the fitted g^ functions and the measured data. The 
deviation A</ 2 )(0) between the fitted value of g^ 2 \0) and the measured datapoints 
is below 0.07 for emitter CI and below 0.15 for emitter C3 without any background 
correction and is mainly governed by the noise of the g^ function. For emitter C4, 
Ag( 2 )(0) = 0.04 is obtained after background substraction. 

Treating the population dynamics of single SiV centres in the framework of a three 
level system marks a significant simplification as we observe four radiative transitions 
at low temperature. Our </ 2 ) measurements correlate photons from all four different 
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Figure 8. Excitation power dependent g' 2 ' functions for emitter CI (a) at room 
temperature and (b) at 30 K. he upper part of each graph shows a large range 
of delay times r to enable observation of the bunching; the lower part displays a 
reduced range of r to reveal the antibunching in detail. 



transitions due to broadband (20 nm) spectral filtering. Thus, we measure an 'effective' 
g( 2 ^ function with properties dominated by the fine structure lines with high relative 
intensity. The full conformity of the measured g^ data with the fitted three level 
g( 2 ) function justifies the description as a three level system and excludes, e.g., that 
several shelving states with different lifetimes have to be considered which would 
induce additional time constants [55]. However, we note that detailed consideration 
of the line intensities in the temperature dependent fluorescence spectra (method see 
[14] ) furthermore indicates that population transfer takes place in the excited state 
(thermalization). Such processes are also not taken into account by the simplified 
three level model. 

From the fits of the </ 2 ' functions, we obtain the power dependent parameters 
a, Ti and t-i in ([3]) summarized in figure [9] The parameter range of a, t± 
and T2 is in accordance with previous room temperature observations [6j 12 lj : 
the comparison of room temperature and low temperature data indicates mostly 
temperature independent population dynamics. The rate coefficients fey (transition 
from state i to j) fully describe the internal population dynamics of the SiV centres. 
In previous work [6, 21 , we developed a model for the internal dynamics: fc2i and /c23 
are intensity independent, fci2 depends linearly on the excitation power P (fei2 = o~P) 
and the de-shelving rate coefficient £31 is given as 
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Figure 9. parameters a, t\ and T2 for (a) emitter CI, (b) emitter C3, (c) 

emitter C4, at room temperature and cryogenic temperature (30K C1,C4; 20 K, 
C3). Filled squares: data points at room temperature (Fit red curve), filled stars: 
cryogenic temperature (Fit blue, dashed curve). For discussion of the employed 
population dynamics model see text. 
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fc 21 (MHz) 


fc 23 (MHz) 


fc^(MHz) 


d (MHz) 


cr(MHz//iW) 


c (/xW) 


RT 
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17.4 
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11.9 


5.21 
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30 K 
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2.79 
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Em C3 


fc 2 i(MHz) 


fc 23 (MHz) 


fe^(MHz) 


d (MHz) 


a (MHz///W) 


c (/zW) 


RT 


770.1 


11.1 


0.749 


5.65 


3.79 


217 


20 K 


889.5 


5.73 
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3.3 
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fc 2 i(MHz) 


fc 23 (MHz) 


fc^(MHz) 


d (MHz) 


a (MHz/VW) 


c (nW) 
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1053.6 
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3.44 


1.83 


201 
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874.8 


18.8 


0.24 


3.4 


11.0 


8.1 



Table 2. Rate coefficients for emitter CI, C3 and C4 deduced from the limiting 
values of a, n and ti using the three level model including intensity dependent 
dc-shclving. 



Equation (01 introduces the possibility that the colour centre is transferred from the 
shelving state (state 3) to the ground state (state 1) as a consequence of laser excitation 
(intensity dependent de-shelving), in addition to the spontaneous relaxation to state 1 
represented by fc^. Figure |H] summarizes the power dependent curves for a, t\ and r 2 
obtained in this model. As discernible from figure |9l the model suitably describes the 
power dependence of a, t\ and r 2 with the exception of r 2 for emitter C3. Thus, an 
intensity dependent de-shelving process is present at room and cryogenic temperature. 

To estimate the rate coefficients in this model, we obtain the limiting values 
t 2 °, t%° and a°° for vanishing and infinite excitation power from the data given in 
figure [5J For details on the model and the procedure to determine the rate coefficients 
see [6l 121). The resulting rate coefficients are summarized in table [2] fc 2 i is mainly 
determined by the limiting value r° which is an estimate for the excited state lifetime 
of the colour centre. As the deviation of low temperature and room temperature values 
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of fc 2 i is smaller than 20 % without a clear trend, we conclude that no temperature 
dependent activation or de-activation of non-radiative channels occurs in accordance 
with the observed brightness of the emitters. 

The emitters reveal significantly different shelving dynamics represented by the 
rate coefficients ^23, k^ and d. However, the room and low temperature parameters 
for the same emitter are in the same order of magnitude, thus the shelving dynamics 
do not fundamentally change upon cooling. For emitter Cf, k^ is inconsistent at 
high and low temperature. However, we note that this might be connected to the 
challenging determination of r° which defines k^: at low temperature, the lowest 
excitation power used was 0.1P sat ; at room temperature 0.02P sat could be employed 
due to the increased performance of the confocal setup at room temperature, and 
thus the reduced measurement time. The observation of a temperature independent 
shelving dynamics excludes thermally activated shelving or de-shelving processes. This 
finding is in contrast to observations on NV centres for which the shelving state is found 
only 37 meV below the excited state and the centres can be thermally re-excited from 
the metastable shelving state at elevated temperature {53] , 

4. Conclusion 

We use cryogenic temperature luminescence spectroscopy to investigate single SiV 
centres in NDs on Ir as well as an ensemble of SiV centres in a high quality, low stress 
homoepitaxial CVD diamond film. As a result of stress in the NDs, the observed 
individual fine structure splittings for single emitters vary significantly. However, we 
find no discernible trend of the modifications. We measure the temperature dependent 
homogeneous broadening which is best described by a T 3 dependence and might be 
associated to phonon-broadening in the presence of impurities. Spectral diffusion 
additionally broadens the emission lines to 25-f60 GHz at around 5 K for most of 
the single emitters. In contrast, the ensemble in the homoepitaxial CVD film shows 
a linewidth of individual fine structure components of down to 9 GHz. The line 
position of the ZPL is found to blue shift by 20-30 cm" 1 upon cooling from room 
temperature to 5 K. The temperature dependence of the shift is consistent with 
a aT 2 + bT 4 dependence and is related to quadratic electron-phonon coupling and 
lattice contraction. However, due to the undetermined stress response of SiV centres, 
we cannot infer the dominant mechanism. All emitters reveal population dynamics 
that are reasonably described in the framework of a three level intensity dependent de- 
shelving model at cryogenic as well as room temperature. The population dynamics 
is mostly temperature independent without de-activation of non-radiative processes 
or changes in shelving or de-shelving dynamics. 
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